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The strange case of a sub-DLA with very little HI 
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Abstract. We report a deep search for HI 21cm emission from the z = 0.00632 sub-DLA toward PG1216+069 with the 
Giant Metrewave Radio Telescope. No emission was detected and our 5cr upper limit on the mass of any associated galaxy is 
Mm 5 10' Mq, nearly 3 orders of magnitude less than M^j. The z ~ 0.006 absorber is thus the most extreme known deviation 
from the standard paradigm in which high column density quasar absorption lines arise in the disks of gas-rich galaxies. 
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1. Introduction 

Neutral atomic gas clouds at high redshift can currently be 
detected only via the absorption lines they produce if they 
happen to lie along the line of sight to an even more dis- 
tant quasar Although present technology allows one to de- 
tect Lyman-a absorption lines from clouds with as low an HI 
column density as Nhi ~ lO'"* cm"^, the bulk of the neutral 
gas at z ~ 3 turns out to be in extremely rare systems, with 
column densities higher by several orders of magnitude, viz. 
Nhi ^ 10^" cm"^. At such large Nhi, the optical depth is sub- 
stantial even in the Lorentzian wings of the Lyman-a profile; 
these systems are called damped Lyman-a absorbers or DLAs. 

In the local universe, the disks of spiral galaxies have 
characteristic HI column densities > 10^" cm"^ over a large 
spatial scale. High z DLAs are thus natural candidates for 
the precursors of today's large galaxies (e.g. I Wolfe 1988t . 
The nature of high redshift DLAs and their evolution with 
redshift have hence been of much interest to astronomers 
studying galaxy formation and evolution. Unfortunately, cos- 
mological dimming makes it impossible to directly image 
high redshift DLAs; the information obtainable is hence 
limited to the gas that lies along the naiTow pencil beam 
that is illuminated by the background quasar It is primar- 
ily for this reason that, despite decades of study, the na- 
ture of high redshift damped absorbers remains an unsettled 
issue, with models ranging from large, rotating disks (e.g. 



IProc haska & Wolfe 1997) to small, merging proto-galaxies 
(e.g. ,Haehnelt, Steinmetz & Ranch 1998^ Detailed absorption 
studies with 10-metre class optical telescopes have established 
that current samples of high redshift DLAs typically have low 
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metal abundances (< 0.1 solar; IPettini et al. 1997t and that 
these abundances evolve slowly, if at all, with redshift (e.g. 
Prochaska et al. 2003 Kulkarni et al. 2004 1. This is somewhat 
surprising if the absorbers are indeed typically the precursors 
of large galaxies like the Milky Way. 

Low redshift (z < 0.5) DLAs offer the possibility of direct 
identification of the absorbing galaxy by optical imaging, (e.g. 
le Brun et al. 1997 Rao et al. 2003 Chen & Lanzetta 2003}- 
Since the HI mass function is dominated by bright galax- 
ies (e.g. I Zwaan et al. 1997t . one might a priori expect that 
low z DLAs should primarily be associated with large spirals. 
Interestingly enough, however, it has been found that a wide 
variety of galaxy types are responsible for damped absorption 
at low redshifts (e.g. Rao et al. 20031 IBowen et al. 20011 . with 
no particular type dominating the cuiTent (admittedly rather 
small) low z sample. One possible explanation is that, although 
the total HI mass is dominated by large spiral galaxies, the HI 
cross-section at the ~ 10^° cm"^ threshold varies across galaxy 
type in such a way as to yield an even distribution of absorbers 
across a range of galaxy luminosities (Zwaan et al. 2002 1. On 
the other hand, it is also possible that, although the optical 
counterparts to these low redshift DLAs are faint, they have 
unusually large HI envelopes and thus, a large total mass. Such 
a preferential selection might well occur since absorption sur- 
veys are biased toward objects with a large HI cross-section. 
Unfortunately, it is very difficult to observationally test the lat- 
ter scenario as this requires detection of the HI 21 -cm line and 
current radio instrumentation limits searches for 21 -cm emis- 
sion from galaxy-sized objects to the very nearby universe, 
z < 0.1. Almost no DLAs are known at these low redshifts 
and hence, searches for 21 -cm emission have been carried out 
in only three absorbers; in all cases, the HI mass has been 
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found to be lower than M*p the typical mass of local bright 
spirals JLane 20001 Ikanekar et al. 20011 IBowen et al. 20011 
[Chengalur & Kanekar"2002f Identifying such z < 0.1 DLAs 
is of much importance as it is only these systems that can be 
followed up in detail in all wavebands, to estimate the total 
mass (as opposed to the luminous mass). We note, in passing, 
that such DLAs cannot be identified (except by accident) us- 
ing the Mgll selection criterion of Rao & Turnshek (2000), as 
the Mgll lines are redshifted into optical wavebands only for 
redshifts z > 0.1. 

Recently, Tripp et al. (2004) reported the discovery of an 
extremely low redshift (z = 0.00632) absorption system to- 
ward the quasar PG1216-I-069, with a high HI column density, 
Nhi~ 2 X 10'*^ cm"^. While the absorption profile shows clear 
damping wings, its column density is somewhat smaller than 
the cut-ofF of 2 X 10^** cm"- that has traditionally been used 
to define DLAs (e.g. lWolfe et al. 1986> . This system is hence 
classified as a sub-DLA. However, it should be emphasized that 
the column density threshold used to define DLAs is entirely 
arbitrary and that the absorber properties do not show any sharp 
qualitative differences here. 

Interestingly enough, although the new absorber is at a very 
low redshift, it has an extremely low metallicity, ~ 1/40 Solar 
( |Tripp et al. 2004} , one of the lowest measured in the gas phase 
in the nearby universe and, in fact, similar to that of high z 
DLAs (e.g.l^ochaska et al. 2003j. No optical counterpart can 
be seen in the Hubble Space Telescope (HST) image of the 
field ( |Tripp et al. 2004 i; the nearest L, galaxy has a projected 
separation of ~ 252 /z^j kpc'. We report here the results of 
a search for HI 21 -cm emission from the z ~ 0.006 sub- 
DLA, using the Giant Metrewave Radio Telescope (GMRT; 
ISwarupetal. 1991> . 

2. Observations, data reduction and results 

The GMRT observations were conducted on the 28th and 30th 
of August 2004. The total on-source time was ~ 8 hours, with 
29 antennas. The observing bandwidth of 2 MHz, centred at 
1411.5 MHz, was divided into 128 spectral channels, yielding 
a spectral resolution of ~ 15.6 kHz. The total velocity coverage 
was thus ~ 425 km s"' , with a resolution of ~ 3.3 km s"' . Flux 
and bandpass calibration were done using short scans on 3C147 
and 3C286 at the start and end of each observing run, while the 
compact source 1 150-003 was used for phase calibration. The 
flux density of 1 150-003 was measured to be 2.7 + 0.1 Jy. 

Data reduction was carried out using standard tasks in 
"classic" AIPS. Visibilities from the different days were cali- 
brated separately and then combined together using the AIPS 
task DBCON; the combined data were used to make the final 
images, with the task IMAGR. We note that the GMRT has 
a hybrid configuration, with fourteen antennas in a central ar- 
ray (the "central square") and the remaining sixteen distributed 
in the three arms of a "Y" ( |Swarup et al. 1991} . The central 
square antennas yield baselines of < 1 km (i.e. U-V coverage 
out to ~ 5kA at 1411.5 MHz), while longer baselines (up to a 

' We assume fi,„ = 0.27, n^ = 0.73 and Hq = 71 km s"' Mpc"' 
throughout this paper. 



Table 1. RMS noise levels and 5cr HI mass limits 



Shpbw RMS rnpBw Mhi 

bmaj"x bmin" mJy/Bm h^l kpc 10* Mj 



39.4 X 37.7 


0.62 


5.0 


11 


26.1 X 23.1 


0.52 


3.3 


8.8 


12.3x9.9 


0.45 


1.6 


7.7 


2.9 X 2.5 


0.30 


0.4 


5.0 



maximum length of ~ 25 km, or a U-V coverage to ~ l20kA at 
1411.5 MHz) are obtained with the arm antennas. Continuum 
images and spectral cubes were hence made with a variety of 
U-V ranges and tapers, allowing a search for HI and contin- 
uum emission at various spatial resolutions ranging from ~ 3" 
to ~ 40". The continuum emission was subtracted using two 
independent approaches: (1) in the U-V plane, using the task 
UVSUB and (2) in the image plane, using the task IMLIN. The 
two strategies gave very similar results; the numbers quoted 
below are for the U-V plane approach. 

Finally, the GMRT does not do on-line Doppler tracking, 
implying that any required Doppler shifts must be applied to 
the data off-line. In the present case, however, it was not neces- 
sary to apply a differential Doppler shift to the data of the two 
observing runs since this shift was found to be small, relative 
to the channel separation. No off-line Doppler corrections were 
hence applied while making the image cubes. 

The central ninety channels of the band were averaged to- 
gether to map the continuum, at various spatial resolutions. The 
low resolution (~ 39.4" x 37.7") image shows good agree- 
ment with the NRAO VLA Sky Survey map of the same region 
JCondonetal. 199^ . Our somewhat higher sensitivity allows 
us to detect faint emission from the quasar itself; we measure 
a flux density of ~ 2.8 + 0.3 mJy. The peak flux density in the 
image is ~ 14.6 ± 0.3 mJy. 

The data cubes were examined for line emission at a variety 
of spectral resolutions; in all cases, no significant line emission 
was found in the vicinity of the quasar line-of-sight. Besides a 
visual inspection, the AIPS task SERCH was used to search for 
line emission in the different cubes. No statistically significant 
emission features were detected in the cube, except from the 
galaxy VCC 297; this system is discussed later. The highest 
resolution (2.9" x 2.5") cube did show a weak (3.3cr) emis- 
sion feature at the quasar location (albeit ~ 70 km s"' away 
from the optical redshift); if real, it corresponds to an HI mass 
of 4 X IQ^Mq. However, this is not very believable as the fea- 
ture is spatially unresolved and seen only in a single velocity 
channel, after smoothing the data to a resolution of 13 km s"^. 
Further, if real, this concentration of gas would imply an HI 
column density > 10^' cm"^, (i.e. considerably larger than the 
column density of the sub-DLA) and it would be surprising that 
no associated metal lines were found at this velocity. 

Tabled summarizes our results, for a representative selec- 
tion of spatial resolutions; the data have been smoothed to a 
velocity resolution of 20 km s"' in all cases, a fairly typical 
velocity width for a small galaxy. The different columns in the 
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table are: (1) the half-power beam width (HPBW) of the syn- 
thesized beam (i.e. the spatial resolution, in arcsecs.), (2) the 
RMS noise (in mJy/Bm) at this spatial resolution, (3) the phys- 
ical distance at z = 0.00632 (in kpc), corresponding to the syn- 
thesized HPBW, and (4) the 5cr upper limit on the HI mass 
of the z ~ 0.006 sub-DLA at this spatial resolution (in units 
of IQ^Mq), assuming that the HI profile has a top-hat shape, 
with a velocity width of 20 km s"'. Note that lower spatial res- 
olutions are obtained by weighting down data from the more 
distant antennas; this implies that the mass limit improves at 
higher spatial resolution. 

3. Discussion 

The 5cr upper limits on the HI mass of the galaxy listed in the 
last column of Table [2 are extremely small, nearly 3 orders of 
magnitude lower than M*j. This is quite remarkable, given the 
high HI column density estimated from the Lyman-a line. Of 
course, it should be pointed out that the mass limits quoted in 
the table are based on the assumption that the absorber is en- 
tirely contained within the synthesized beam. It is, in principle, 
possible that the HI emission is actually spread over a larger 
angular scale and is not detected in our interferometric observa- 
tions either because the flux per synthesized beam falls below 
our detection threshold or because the flux is resolved out or 
due to a combination of both effects. We initially consider this 
possibility, before discussing the implications of our results for 
the nature of the z = 0.00632 sub-DLA. 

If the absorbing galaxy were far larger than the GMRT 
synthesized beam, the mass limits of Table ^ would re- 
fer only to the HI mass within the beam and not to that 
of the entire galaxy. For example, if the galaxy were a 
face-on disk of diameter 10 kpc, twice our coarsest res- 
olution, the upper limit to its total HI mass would be a 
factor of ~ 4 larger than that listed in the first line of 
Table ^ i-S- Mm < 4 x 10^ Mq. However, dwarf irregular 
galaxies with low HI masses (< lO^M©) have fairly small 
HI diameters (e.g . Stil & Israel 2002a Stil & Is raer2002bl 
[Begum, Chengalur & Hopp 2003|) . For example, all ten dwarf 
irregulars with Mhi < 10** Mq in the sample of Stil & Israel 
(2002a, 2002b) have size < 5 kpc, our coarsest resolution. In 
fact, of the 27 dwarfs in the Stil & Israel sample with estimates 
of HI extent, the only galaxies that have physical sizes larger 
than 11 kpc are the five systems with Mm > 6.5 x 10^ Mq (i.e. 
considerably larger than the upper limit in Table^J. Further, if 
we assume that Nhi ~ 2 x 10'^ cm"^ throughout the absorber, 
the total HI mass in an area 5 kpc in diameter is ~ 3 x 10'' M©, 
comparable to our detection threshold. The latter is a severe 
under-estimate of the true HI mass, if the absorption arises in a 
normal galaxy, since large fractions of even dwarf galaxy disks 
have Nhi <: 10^'^ cm"^. It is thus highly unlikely that the absorp- 
tion arises in a dwarf galaxy whose HI disk is so extended that 
the flux per synthesized beam falls below our detection thresh- 
old. Similarly, the second scenario, that the flux is resolved out 
because it arises from a highly uniform HI distribution, is also 
rather improbable as this would require a very small velocity 
gradient over the entire spatial extent of the gas, for emission 
to not be detected even at the highest velocity resolutions. 



The nearest known Hl-rich galaxy to the sub-DLA is 
VCC 297, which has L = 0.25L* ( |Impey et al. 1999) , 
Mhi - 2.7 x IO^Mq and a velocity width W50 of 
145 km s ' ( |Giovanelli, Avera & Karachentsev 1997 ;); this is 
~ 87 h'l kpc away from flie QSO sighthne. VCC 297 lies at 
the edge of our field of view, where imaging is particularly dif- 
ficult because of the asymmetry of the edges of the GMRT pri- 
mary beam. The emission velocity of VCC 297 also lies at the 
edge of the observing band where the sensitivity is lower; fur- 
ther, there appears to be some low level interference at these 
edge channels (note that none of these problems are present 
at the centre of our field of view and observing band, where 
the emission from the sub-DLA would arise). Next, the optical 
extent of VCC 297 is a factor of ~ 3 larger than our synthe- 
sized beam; its HI extent is likely to be at least this large, re- 
ducing the flux per synthesized beam. Despite these problems, 
HI emission from VCC 297 is reliably detected in the spec- 
tral cubes (see Fig. The emission was readily apparent to 
the eye and was also detected by the automatic search algo- 
rithm SERCH, at a signal-to-noise ratio of ~ 12. Finally, as a 
concrete counter-example, a dwarf galaxy at z = 0.005 1 with 
Mhi ~ 2 X 10^ Mq was serendipitously discovered in a diff'er- 
ent GMRT observation, with essentially identical observational 
setup and RMS noise ( [Chengalur et al. 2004 1. From all of the 
above, we conclude that it is highly unlikely that the sub-DLA 
arises in a galaxy with HI mass substantially larger than the 
limits quoted in Table [2 Further, if one assumes a uniform HI 
columndensity of 2x10''' cm"- throughout the absorber (and a 
disk geometry), our HI mass limit of ~ IQ^Mq implies an upper 
limit of ~ 9 kpc to the size of any associated galaxy. It is still 
not impossible, of course, that the sub-DLA arises in a more 
massive, highly extended, smooth HI cloud, but this would re- 
quire it to be completely unlike any known low redshift galaxy. 

As discussed in detail by Tripp et al. (2004), the z ~ 0.006 
sub-DLA is highly unusual in having a very low metallicity 
([O/H] - -1.6;^Q jj) despite being at low redshift and lying 
in a region of high local galaxy density (the outskirts of the 
Virgo cluster). The sub-DLA also shows an under-abundance 
of nitrogen and an over-abundance of iron (implying a lack of 
dust) which also argues in favour of a system that is at an early 
stage of chemical evolution (Tri pp et al. 2004| . These abun- 
dance patterns are very different from those characteristic of 
local spirals, and are consistent with the absence of a lumi- 
nous optical counterpart to the sub-DLA in the HST image. Of 
course, the lack of a detected optical counterpart does not en- 
tirely exclude a scenario in which the sub-DLA arises in a mas- 
sive galaxy, since such a system might well be hidden beneath 
the point spread function of either the QSO itself or a nearby 
foreground star located near the quasar (Tripp et al. 2004 1; al- 
ternatively, as discussed in the introduction, the absorber might 
be anomalously gas-rich for its luminosity. However, the strin- 
gent upper limit on the HI mass of the absorber provided by 
our observations indeed rules out this possibility. 

The abundance patterns and low metallicity of the z ~ 
0.006 absorber are consistent with an origin in (i) a low metal- 
licity, gas-rich, blue compact dwarf (BCD) like I Zw 18 or 
SBS0335-052, (ii) a high velocity cloud (HVC), (iii) a dwarf 
spheroidal galaxy or (iv) a dark mini-halo left over from the 
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Fig. 1. Integrated HI 21 -cm emission profile toward VCC 297, 
made from the 39" x 34" resolution image cube. The contours 
are at 0.08,0. 12,0. 17,0.25 and 0.36 Jy/Bm km s"' . 



epoch of reionization (see Tripp et al. (2004) and references 
therein). However, the mass limits in Table [2 are significantly 
lower than HI masses typical of low metallicity BCDs (> 
10^Mo;[van Zee et al. 1998 Pustilnik et al. 2001 ). It is also un- 
likely that the sub-DLA is associated with an HVC, given that 
the nearest L* galaxy (NGC4260) is ~ 252 /i^j' kpc away (in 
projection) and HVCs have not been detected at such large 
distances from their parent galaxies, despite deep searches 
(e.g. Pisano et al. 2004 1. In fact, for M31 (the only galaxy that 
has been sensitively searched for large-scale diffuse HI emis- 
sion), the typical HI column density falls to values smaller 
than lO'^ cm~^ by the time one reaches galacto-centric dis- 
tances of 200 kpc (Braun & Thilker 2004). Similarly, Miller & 
Bregman (2004) place a limit of ~ 25 kpc on the distance of 
HVCs from their parent galaxy; this implies that the absorp- 
tion is also unlikely to arise from an HVC associated with the 
fainter galaxy VCC 297, which is ~ 87 h^^ kpc away. The most 
likely counterpart is hence an extreme dwarf galaxy (or a dark 
mini-halo), which is surprising, given the relatively small size 
of these systems. How likely is it that such a system would 
show up in an unbiased search for absorption? While there 
are no existing measurements of the typical size of faint dwarf 
galaxies at HI column densities of 2 x lO'^ cm"^, a linear ex- 
trapolation of the results of Zwaan et al. (2002), for systems 
with higher mass and column density, (and further assuming 
the faint end slope of the HI mass function to be -1.2) indi- 
cates that the HI cross-section offered by ~ 10^ M© galaxies 
is significantly more than an order of magnitude smaller than 
that of M*j galaxies. The implied low probability of a sub-DLA 
arising in the former class of systems may suggest that our un- 



derstanding of the HI mass function is incomplete at the low 
mass end. 

In summary, contrary to the a priori expectation that sys- 
tems found via absorption line searches should be biased to- 
ward galaxies with large HI mass, our upper limit to the to- 
tal HI content of the z ~ 0.006 sub-DLA is ~ three orders of 
magnitude lower than M*j. This is by far the smallest HI mass 
Umit placed on gas associated with DLAs or sub-DLAs. The 
z = 0.00632 sub-DLA toward PG1216+069 thus appears to be 
the most extreme known deviation from the standard paradigm 
for high column density absorbers. 
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